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I 

The  development  of  the  DOPPLR  program  was  initiated  in  1969  to  support  the  U.S. 
Army  Engineer  Topographic  Laboratories  in  its  analysis  of  the  Doppler  Backpack 
Positioning  Equipment  with  the  knowledge  that  the  Gcoceiver  was  being  designed  as 
the  prime  Department  of  Defense  geodetic  satellite  tracking  system.  The  program  has 
^ evolved  into  the  main  Defense  Mapping  Agency  Topographic  Center  (DMATC)  point 

. positioning  routine,  and  has  been  used  extensively  not  only  for  production  work,  but 

also  in  conducting  numerous  equipment  tests  and  in  the  analysis  of  the  Navy 
Navigation  Satellite  System. 

The  willing  and  capable  assistance  of  DMATC  employees  in  the  Department  of 
Geodesy,  too  numerous  to  mention,  is  hereby  acknowledged. 
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DOPPLR-A  POINT  POSITIONING  PROGRAM 
USING  INTEGRATED 
IX)PPLER  SATELLITE  OBSERVATIONS 


SECTION  I.  INTRODUCTION 

1.  DATA  REDUCTION.  </.  This  report  is  concerned  with  the  point  positioning  mode 
of  reduction  of  Doppler  satellite  tracking  data.  It  describes  the  mathematical 
development  of  the  pertinent  equations  and  describes  the  operation  of  a computer 
program,  DOPPLR,  for  accomplisliing  the  data  reduction. 

b.  Point  positioning  uses  the  Doppler  shift  of  stable  signals  broadcast  from  satellites 
to  determine  the  precise  geodetic  coordinates  of  the  receiving  antenna.  Thus,  these 
coordinates  are  the  principal  unknowns  of  the  problem,  and  all  known  error  sources 
are  accounted  for  in  the  mathematical  model  in  order  to  obtain  the  greatest  possible 
accuracy  and  precision.  Two  features  distinguish  point  positioning  from  other  data 
reduction  methods.  First,  the  satellite  ephemeris  is  obtained  from  some  outside  source 
and  is  held  fixed  in  all  computations:  second,  the  solution  is  made  in  a multipass  mode, 
involving  the  simultaneous  adjustment  of  data  from  many  passes  of  one  or  more 
satellites  over  a single  receiving  station. 

2.  DOPPLR  PROGRAM,  a.  As  a diagnostic  aid,  the  DOPPLR  program  will  perform 
navigation  mode  solutions,  solving  for  the  receiving  station' coordinates  separately  for 
each  satellite  pass.  The  program  may  also  be  used  in  an  error  analysis  mode.  In  this 
mode,  simulated  data  are  generated  and  perturbed  by  specified  random  and  bias  errors. 
The  perturbed  data  are  then  adjusted  to  show  the  effects  of  the  specified  error  sources. 

b.  The  program  operates  on  integrated  range  rate  data  produced  by  Geoceiver,  ITT 
5500,  and  Backpack  observing  equipment.*  Each  data  point  consists  of  an  observed 
Doppler  count,  the  epoch  of  the  beginning  of  the  count,  and  the  time  interval  over 
which  the  count  was  accumulated.  The  program  applies  necessary  corrections  for 
refraction  and  other  effects.  A separate  preprocessor  may  be  used  to  put  data  recorded 
by  various  receivers  into  the  proper  format  for  operation  of  the  program. 

SECTION  II.  MATHEMATICAL  DEVELOPMENT 

3.  BASIC  CONCEPTS,  u.  Most  Doppler  receivers  work  in  the  same  basic  way:  the 
received  frequency,  /^,  is  mixed  with  a locally  generated  reference  frequency,  jp,  to 
produce  a beat  frequency,  ff^  . The  cycles  of  this  beat  frequency  arc  counted  from 


* Any  mention  herein  of  a commercial  product  does  not  constitute  endorsement  by  the  U.S.  Government. 
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time  7 1 to  time  Tj  . The  Doppler  count  (often  called  integrated  Doppler  count)  is 
thus: 


The  interval.  = Tj  — 7'i , over  which  the  Doppler  count  is  accumulated,  depends 
on  the  type  of  receiver  being  used  and  the  mechanization  used  to  produce  the  time 
marks  7',  and  7'j  . The  intervals  used  by  existing  receivers  range  from  less  than  1 
second  to  approximately  2 minutes.  The  average  beat  frequency  over  the  interval  is 
N/A7'.  When  the  interval  is  sufficiently  short,  this  quantity  may  be  regarded  as  an 
approximation  of  the  instantaneous  beat  frequency  and  may  be  ultimately  related  to 
range  rate. 

b.  In  the  EX)PPLR  program,  the  observation  equation  is  always  written  in  terms  of 
the  Doppler  counts,  which  can  be  ultimately  related  to  range  difference.  The  exact 
form  of  the  observation  equation  depends  on  the  receiver  and  the  source  of  the  time 
signals.  Two  sources  of  time  marks  are  recognized  by  the  program. 

(1)  The  satellite  clock  is  the  source  of  the  time  signals  when  Navy  Navigation 

Satellites  (NNS  |or  Navsats] ) are  tracked.  These  satellites  transmit  stable  carriers  at 
approximately  400  and  150  megahertz  (MHz).  Digital  data  are  phase  modulated  on  the 
carrier  at  the  rate  of  6,103  bits  per  2 minutes  of  satellite  clock  time.  The  satellite  clock 
is  generally  kept  within  50  microseconds  (pis)  of  Universal  Time  Coordinated  (UTC)  by 
daily  monitoring.  The  format  of  the  digital  data  repeats  every  2 minutes,  and  the 
beginning  of  the  format  (satellite  time  mark)  may  be  recognized  by  a unique  bit 
pattern.  Thus,  2-minute  intervals  may  be  realized  simply  by  recognizing  successive 
satellite  time  marks.  Shorter  intervals  may  be  realized  by  counting  the  bits  within  the 
2-minute  format.  When  the  appropriate  number  of  bits  have  been  counted  by  the 
receiver,  a time  mark  is  sent  to  the  Doppler  counting  circuits.  In  the  receivers 
considered  here,  the  accumulation  of  the  current  Doppler  count  ends  and  a new  one  . 

begins  at  the  first  positive-going  zero  crossing  of  the  beat  signal  following  the  reception 

of  the  time  mark.  Thus,  no  partial  cycles  of  the  beat  signal  are  lost  between  successive 

counts,  and  the  counts  are  called  continuously  counted  integrated  Doppler.  < 

(2)  The  Doppler  beacons  in  the  GEOS  (Geodetic  Earth-Orbiting  Satellite)  scries 
of  satellites  transmit  pure  tones  at  approximately  324  and  162  MHz.  In  this  case,  the 
time  marks  are  generated  by  the  local  clock.  As  before,  the  Doppler  counts  begin  and 
end  at  the  first  positive -going  zero  crossing  of  the  beat  signal  following  the  reception  of 
the  locally  generated  time  mark. 
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c.  For  both  the  NNS  and  the  GEOS  satellites,  the  measurement  of  the  Doppler 
count  is  made  on  the  higher  frequency,  while  the  lower  frequency  is  used  only  for  the 
computation  of  ionospheric  refraction.  Both  frequencies  of  a pair  are  derived  from  the 
same  basic  oscillator  and  are  coherent  at  transmission.  The  actual  transmitted 
frequencies  are  offset  from  the  nominal  values  for  both  satellite  types.  The  offset  for 
the  NNS  is  —80  parts  per  million  (ppm),  or  —32  kilohert;^  (kHz)  at  400  MHz;  for  the 
GEOS  satellites,  it  is  —50  ppm,  or  —16.2  kHz  at  324  MHz. 

d.  For  the  receivers  considered  here,  the  locally  generated  reference  frequencies  are 
the  nominal  frequencies  of  400  MHz  for  the  navigation  satellites  and  324  MHz  for  the 
GEOS  satellites.  The  Doppler  shift  for  an  observer  fixed  on  the  Earth  never  exceeds 

±20  ppm  of  the  transmitted  frequency.  Thus,  the  received  frequency  (transmitted 
frequency  plus  Doppler  shift)  is  always  lower  than  the  reference  frequency,  and  the 
beat  frequency  is  always  in  the  sense  jy  =/q  — /^  . 

4.  BASIC  FORMULAS 

d.  Sdvij^ation  Satellite  Keccptioti 

(1)  Counting  Intervals.  For  the  receivers  considered  here,  the  first  time  mark  is 
generated  upon  recognition  of  the  first  format  beginning  (satellite  time  mark) 
following  signal  acquisition  and  receiver  lock.  Successive  time  marks  are  generated  as 
follows. 


(a)  Backpack— A time  mark  is  generated  every  6,103  bits,  or  2 minutes  of 
satellite  time. 

{bj  ITT  5500— Each  2-minute  interval  is  divided  int<j  26  subintervals.  The  first 
25  subintervals  are  each  234  bits,  or  4.601016  seconds,  of  satellite  time  in  length.  The 
26th  subinterval  lasts  for  253  bits,  or  4.974603  seconds  of  satellite  time. 

(c)  Geoceiver— Each  2-minute  interval  is  divided  into  4 subintervals  as  follows: 


.Subintcrv;tl 

Number  of  Bits 

Time  Duration 

1 

1638 

32.207111  sec. 

2 

1404 

27.606094  sec. 

3 

1638 

32.207111  sec. 

4 

1423 

27.979683  sec. 

Thus,  in  all  cases,  the  time  and  duration  (at  the  satellite)  of  each  Doppler  count  is  fixed 
and  known  simply  from  the  position  of  the  Doppler  count  in  the  format  output  by  the 
receiver.  Additionally,  the  Geoceiver  outputs  the  reading  of  its  local  clock  at  the 
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beginning  and  the  end  of  each  Doppler  count.  From  these  data,  the  epoch  error  of  tlic 
local  clock  with  respect  to  the  satellite  clock  may  be  computed.  This  feature  also 
allows  the  computation  of  a small  quantity  known  as  the  “partial  cycle  count 
correction.”  This  small  correction,  which  cannot  be  computed  for  data  received  by 
other  receivers,  eventually  accounts  for  the  higher  precision  and  accuracy  of  the 
Geoceiver. 

(2)  The  Doppler  Equation,  {a)  Let  t be  the  time  of  transmission  from  the 
satellite  of  a time  mark  (or  a bit  which  will  generate  a time  mark  in  the  receiver).  Then 
the  counting  of  the  beat  signal  begins  at  f + At  + e,  where  At  is  the  propagation 
delay  of  the  signal  and  c is  a further  delay  composed  of  the  time  required  to  sense  the 
time  mark  in  the  receiver  and  the  wait  for  the  next  zero  crossing  of  the  beat  frequency. 
The  Doppler  equation  is  written  as 

ti  + Atj  + Cj 

N = r f^dt  . 

ti  + At,  + Cl 


Ignoring  refraction  effects  (which  are  treated  separately),  At,  and  At2  are  given  by 
r, /c  and  rj/c,  where  r,  and  T2  are  the  ranges  to  the  satellite  at  times  t,  and  tj , 
respectively,  and  c is  the  vacuum  speed  of  light.  The  ranges  are  in  turn  given  by  r,  = 
x/(K,  — x)~  and  rj  = V (Rj  — x)^  where  Ri  and  Rj  are  the  (assumed  known) 
position  vectors  of  the  satellite  (in  an  Earth-fixed  basis)  at  t,  and  tj , respectively, 
and  X is  the  receiving  antenna’s  position  vector  (assumed  constant,  but  unknown). 


{b)  Substituting  for  /jj , 

<2  + Af2  +62 

fj  + Afj  + C; 

N = j 4* 

f • 

f,  +Af,  +6, 

t,+  Af,  + 6 

(1) 


Since  /g  is  constant  (or  at  least  can  be  considered  to  be  so  over  the  span  of  a satellite 
pass),  the  first  integral  is  simply  /g(^j  — tt  + A/j  — ■At,  + C2  — e,  ) . The  second 
integral  is  more  complicated,  if  the  delay  e is  ignored,  it  may  be  evaluated  by  arguing 
that  the  number  of  cycles  received  between  f,  +At,  and  f2  +Af2  must  be  the  same 
as  the  number  of  cycles  transmitted  between  f,  and  f2 . or  ff{ti  — /,  ),  where  is 
the  transmitted  frequency.  This  leads  to 


N = foiti  - f,  + Af2  - At,  ) -/,(fj  - G) 

= {fo  -f,)(h  -G)+/o(Ai2  - Af,) 
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or 

N = A/-A7; +/„(A/j  -Af,), 

with 

A/'  = fo  and  A7'^  = fj  - t,  . 

This  form  of  the  Doppler  equation  is  often  used  in  navigation  applications.  A more 
careful  evaluation  of  the  second  integral  in  equation  (1),  considering  the  delay  (e),  is 
gtven  in  paragraph  7.  Using  the  result  given  there, 

^ = /o(*j  - t,  + At,  - At,  +f,  _ r,  +e,  _e,  + i -e,)  + ^(r,  _ r,  )? 

Regrouping  terms,  this  equation  is  rewritten  as 

,\  = A/  A / j + /o  (Afj  - At,  ) + A/-(c2  - f I ) + X («2  - f 1 ) + (tj  — r,  )f  , (2) 

where  AF  = fo  — f(  is  the  frequency  offset  (close  to  32  kHz),  A7'^  = — f,  is  the 

time  interval  at  the  satellite,  and  X = c//^  is  the  wavelength  of  the  transmitted  signal. 

(c)  The  term  (Afj  — Af, ) contains  the  contribution  of  the  Doppler  shift  to 
the  count.  For  convenience,  the  term  is  rewritten  as 

A/-(Af2  - Af,  ) + /^(Afz  - Af,  ) , 
or 

A/'(rj  - r,  )/c  + -^(r2  - r,  ) . 

The  first  term,  A7''(r2  — r,  )/c,  is  labeled  CT  . In  the  program  it  is  computed  from 
approximate  values  of  A/-',  t2  , and  r,  . The  last  term  now  shows  the  Doppler  shift  as 
occurring  on  the  transmitted  rather  than  on  the  reference  frequency.  In  addition,  the 
equation  is  now  in  a form  that  can  be  applied  to  GEOS  satellite  observations. 

{(/)  The  sum  A/-(c2  — Ci ) r(c2  — Ci ) is  labeled  the  PC  (for  partial  cycle) 

correction.  The  PC  term  can  be  identified  as  the  average  beat  frequency  times  the 

difference  in  delays  at  the  beginning  and  end  of  the  counting  interval. 

/<;  = (AF  + if)(e2  -e.)  = fye,  - e,  ) • 

Although  the  individual  delays  are  unknown  and  cannot  be  determined  directly  from 
the  data,  the  difference  62  — e,  can  be  determined  by  comparing  the  counting 
interval  as  measured  by  the  local  clock  and  by  the  satellite  clock.  For  Geoceiver  data, 
this  difference  is  computed  from  approximate  values  of  and  r,  , using 

ej  - e,  = A7'^  - A /'^  - (r,  - r,  )/c  , 
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where  A 7^,  is  the  actual  counting  interval  measured  by  the  (Jeocciver  clock.  The 
average  beat  frequency  for  the  interval  A /'^  is  estimated  by 

- observed  Doppler  count 

fb  = 

For  data  from  other  receivers  not  having  a local  clock,  the  PC  term  is  ignored. 

(e)  The  final  form  of  the  Doppler  equation  is  thus: 

A’  = AFAT^  + Cl  + ^(rj  - r,  ) + I^:  + {(fj  - r,  )e  ■ (3) 


b.  GUOS  Satellite  Reception 

(!)  Counting  Intervals,  (u)  For  GEOS  satellites,  time  marks  are  controlled 
entirely  by  the  local  clocks.  Counting  of  the  beat  frequency  begins  at  the  first  even 
2-minute  time  mark  (according  to  the  local  clock)  following  signal  acquisition  and 
receiver  lock. 

[h)  The  ITT  5500  simulates  the  bit  stream  from  the  satellite  by  the  local 
generation  of  6,103  bits  per  2 minutes  of  local  clock  time.  The  subintervals  within  the 
2-minute  format  are  the  same  as  for  the  navigation  satellites. 

(c)  The  Geoceiver  generates  a time  mark  every  30  seconds,  on  the  minute  and 
the  half-minute,  according  to  the  local  clock.  Again,  the  Geoceiver  also  outputs  the 
reading  of  the  local  clock  at  the  beginning  and  end  of  each  Doppler  count.  These  times 
differ  frotn  even  Vi-minute  times  by  the  wait  for  the  next  positive-going  zero  crossing 
of  the  beat  frequency. 

(2)  The  Doppler  Equation.  Equation  (2)  also  applies  to  GEOS  satellite  reception, 
except  that  the  receiver  delays  do  not  enter  into  consideration.  Since  the  interval  at  the 
satellite  is  not  known,  the  Doppler  equation  must  be  written  in  terms  of  the  interval 
measured  on  the  ground  clock.  This  is  accomplished  by  rearranging  equation  (2)  as 

/V  = A7-A7^+-^(r,  -r,).  (4) 

(3)  General  Form.  The  DOPPLR  program  uses  the  following  single  form  to 
include  equations  (3)  and  (4); 

N = M-^T  + CI  + IK:  + \(r2  - r,  ) + -^  (f:  - G )e  . 

A A 

when  navigation  satellites  are  tracked,  A7'  is  A7'^  and  Cl  is  computed  from 
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approximate  values.  When  GEOS  satellites  are  tracked,  AT  is  and  Cl  = 0 . Tin- 
term  PC  is  zero  in  all  c.tses  except  for  Geoceiver  tracking  of  navigation  satellites.  The 
last  term  is  the  time  bias  unknown  and  is  explained  in  paragraph  7.  When  used  as  the 
basis  for  the  adjustment  of  the  station  location  in  the  point  positioning  mode,  the 
quantities  X,  AT,  R,,  Rj,  R|,  and  Kj  are  always  assumed  perfectly  known.  The 
frequency  offset  AT  is  an  unknown  whose  value  varies  from  pass  to  pass.  The  station 
position  X,  as  contained  in  the  factor  (rj  — r, ),  is  usually  treated  as  completely 
unknown.  However,  provision  is  made  to  constrain  any  combination  of  latitude, 
longitude,  and  height  to  agree  with  input  values;  the  degree  of  the  constraint  is 
determined  by  a priori  input  estimates  of  accuracy.  The  time  bias  unknown  e may  be 
allowed  to  vary  from  pass  to  pass,  satellite  to  satellite,  receiver  to  receiver,  or  it  may  be 
dropped  from  the  equation  and  ignored  altogether. 

5.  C!ORRECTIONS  TO  OBSERVED  DATA.  The  Doppler  count  N of  the  preceding 
section  is  the  ideal  count,  which  would  be  obtained  in  the  absence  of  physical  effects 
on  the  signal.  The  actual  measured  count,  , is  corrected  to  the  ideal  count  N by 

N = d^  — corrections, 

where  the  corrections  are  due  to  tropospheric  and  ionospheric  refraction  and  to  the 
relativity  effect. 


a.  I roposphcric  Refraction.  (1)  The  correction  for  the  tropospheric  refraction 
(TROP)  to  the  Doppler  count  is  computed  using  the  tropospheric  refraction  correction 
for  each  slant  range,  if  the  calculated  range  delay  for  slant  range  t|  is  1 Ri  and  for  rj 
is  77^2,  then 

N + TROP  = A/-A7'  + -^(rj  + 7 Kj  - r,  - TR , ) . 

Q>nsequently, 

TROP  = \(TR^  - TR,). 

I Rj  is  calculated  by  integrating  the  atmospheric  refractivity  along  a line  from  the 
station  to  the  satellite.  The  integral  is  computed  for  both  the  dry  and  wet  components 
of  the  refractivity  and  the  results  added  to  give  the  total  tropospheric  delay  in  meters. 
The  integral  computed  for  each  refractive  component  is 


TR.  = 

I 


R'l'  + upper 

f 

J vV  + RT’eos’e 

RT  + lower 


(5) 
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where 


5 - distance  from  the  center  of  the  Earth  to  a point  on  the  straight  line  propa- 

gation path. 

upper  = upper  limit  of  integration  (height  above  station), 
lower  = lower  limit  of  integration  (height  above  station), 
cos  e = cosine  of  elevation  angle. 

N(6)  = refractivity  at  height  6. 

Rf  = distance  from  Earth’s  center  to  station. 


The  two  quartic  refractivity  model  developed  by  Hopfield'  is  used  for  the  refractivity 
profile 

N(d)  = 10-*  N ' , 


(upper  — lower)^ 


where  N is  the  surface  refractivity. 


(2)  For  Backpack  and  ITT  5500  data,  the  surface  wet  bulb  (T^)  and  dry  bulb 
(Tj)  temperatures  and  pressure  (;>)  are  measured  for  each  pass.  The  surface  refractivity 
is  calculated  from  these  values  using 


77.62  .. 


= [< 


37.19 

7;,/100)^ 


I j Vap  . 


Vap  = Sat  - 0.00066(1.0  0.00115  7’^)  p (Tj  - T ). 


Sat  = 0.4629  X ]0-'«  /-^«  + 0.155  X lO-’  /'^’  + 0.2254  X 10'’ T ’ 
+ 11.9792  X lO^T^’  + 10.7253  X 
+ 33.337  X 10-2  7'^ -K  4.58  . 


(3)  When  Geoceiver  observations  are  processed,  the  surface  refractivity  is  calcu- 
lated from  surface  temperature  (Tj),  pressure  (p),  and  relative  humidity  (77)  values. 

N,  = UAp  . 

' (i 


N, 


= * //  * 0.01  ♦ <'(-37.2465  + 0.213177  y -0.00026rj)j 


The  integral,  equation  (5),  is  evaluated  by  using  the  series  expansion  given  by 
Yionoulis.  ® ^ 

' HOPMKLl)  H.S.  “A  Two.(^iartic  Tropospheric  Rcfraclivity  Profile  for  CorrectinK  .Satellite  Data."  Presented  to 
the  InternaMorial  Symposium  on  Electromagnetic  Distance  Measurement  and  Atmospheric  Refraction,  Boulder 
^ (.olorado  |une  1969. 

^ ^ Tropospheric  Refraction  Effects  on  Range  Measurements."  .4pn/,>,( 

Physu,  laboratory  T,ch„u:al  M,morand„m  TG  1125.  Silver  Spring,  Md.:  Johns  Hopkins  University  July  1970. 
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h.  lonosjilicric  Rcjraction.  in  addition  tt)  tlic  basic  frequency  from  which  the 
Doppler  count  is  measured,  botl)  navigation  and  GEOS  satellites  transmit  on  a second, 
coherently  related  frequency.  A comparison  of  the  reception  of  the  two  frequencies 
yields  a measurentent  of  the  first-order  effects  of  ionospheric  refraction.  Both  the 
Geoceiver  and  the  ITT  5500  produce  a refraction  count  derived  from  the  two 
frequency  observations.  The  ionospheric  correction  (ION)  is  obtained  by  scaling  the 
measured  refraction  count.  The  scaling  factors  for  the  Geoceiver  (Stanscll.  et  al  ^ ) and 
the  ITT  5500.  are 


Frequency  Pair  (Jcocciver  ITT  5500 

400/150  MHz  6/55  24/55 

324/162  MHz  1/9  8/9 

In  the  case  of  the  Backpack,  the  ionospheric  correction  is  applied  to  the  measured 
Doppler  count  by  analog  means  in  the  receiver,  so  that  no  further  correction  for 
ionospheric  refraction  need  be  made. 

r.  Rclatii'ity.  Because  of  special  relativity,  the  frequency  of  the  signal  transmitted 
from  the  satellite  appears  to  an  observer  fixed  on  the  Earth  to  be  lowered  by  a factor 
(1  — rVe^)''’  , where  i’  is  the  velocity  of  the  satellite  in  an  Earth-fixed  frame.  This 
amounts  to  approximately  0.1  Hz  at  the  400  MHz  frequency,  or  3 counts  over  a 
30  second  counting  interval.  This  effect  is  indistinguishable  from  a real  drift  of  the 
satellite  oscillator.  Since  the  program  always  solves  for  the  frequency  offset  (difference 
between  transmitted  and  local  reference  frequency)  on  a pass-by-pass  basis,  it  is  not 
necessary  to  make  a separate  correction  for  the  relativity  effect.  Thus,  the  frequency 
offset  obtained  in  the  solution  is  actually  a composite  of  the  actual  frequency  offset 
and  the  relativity  effect. 

(/.  Corrected  Doppler  liquation.  The  final  form  of  the  observation  equation,  with 
all  known  terms  on  the  left  and  all  unknown  terms  on  the  right,  is 

= d„  - TROr  - ION  - CA  - PC  = A/  A7  + ^(rj  - r,  ) -F (rj  - n )e  . (6) 

6.  LO(!AI,  CLOCK  SYNCHRONIZATION,  a.  The  quartz  crystal  oscillators  used  in 
Doppler  receivers  exhibit  excellent  short  term  stability,  but  typically  drift  slightly 
away  from  their  intended  frequency  over  periods  of  days  or  weeks.  The  receiver  clock 
is  normally  started  at  the  beginning  of  an  operational  period  and  allowed  to  run 
continuously  for  the  whole  operation,  so  that  the  epoch  carried  by  the  local  clock 
typically  drifts  slowly  away  from  the  inital  synchronization  to  UTC. 


STANSELL.  T.A.,  ct  al.  "(Icoceivcr:  An  integrated  Doppler  Geodetic  Receiver,”  .Applied  Physics  Laboratory 
/echtiical  Mt-morandum  TG  710.  Silver  Spring,  Md.:  Johns  Hopkins  University.  July  1965. 


13 


4 


b.  When  navigation  satellites  arc  tracked,  the  local  clock  epoch  error  may  be 
determined  by  comparing  the  recorded  time  of  the  beginning  of  a Doppler  count  with 
the  known  time  (in  UTC)  of  the  transmission  of  the  corresponding  bit  in  the  navigation 
satellite  message,  if  Ar^  is  the  correction  tor  the  local  clock  epoch  error,  the  time 
recorded  by  the  local  clock,  and  t the  corresponding  satellite  (approximately  UTC) 
time,  then  the  UTC  of  the  bcginiiiiig  of  the  Doppler  count  is 


so  that 


f,  + At^  = t + ric  + e , 
A f ^ = t + r/c  + e — . 


With  the  (ieoceiver,  a reading  of  the  local  clock  is  obtained  with  each  Doppler 
count.  The  reading  is  recorded  to  a precision  of  4 /ts,  the  transmission  time  t is 
known  to  within  50  ;js  in  UTC  (the  accuracy  of  the  satellite  clock  in  UTC),  and  the 
propagation  delay  r/c  may  be  predicted  to  a few  ns  . In  the  program,  the  quantity 
Af^  — e is  computed  for  Geoceiver  observations  of  navigation  satellites,  and  an 
average  is  taken  over  all  points  in  the  pass.  Thus,  the  computed  correction  for  the  local 
clock  epoch  error  is 

® ■ i I, 

Assuming  the  partial  cycle  wait  for  the  zero  crossing  of  beat  frequency  to  be  random, 
the  average  value  e of  e may  be  identified  as  the  average  receiver  delay.  The  quantity 
CF.  is  computed  by  the  program. 

r.  In  norma]  field  operations,  the  tracking  of  navigation  satellites  is  the  only  means 
of  determining  the  local  clock  epoch  error  and  calibrating  the  drift  of  the  local  clock. 
Thus,  navigation  satellites  must  be  tracked  for  timing  purposes  even  when  CiEOS 
satellites  are  tracked  for  positioning  or  orbit  determination  purposes.  When  GEOS 
tracking  data  arc  processed  by  the  DOPPLR  program,  the  local  clock  epoch  error  and 
drift  rate  at  a specified  epoch  arc  input  and  used  to  correct  the  recorded  time  of  each 
data  point  to  UTC. 

7.  RECEIVER  DELAY,  a.  The  delay  e is  composed  of  the  wait  for  the  zero  crossing 
of  the  beat  signal  and  the  actual  receiver  delay,  which  is  the  time  required  for  a given 
bit  in  the  satellite  message  to  be  recognized  by  the  receiver  logic  and  a signal  sent  to 
the  Doppler  counting  circuitry.  Let  t be  the  time  of  transmission  of  a bit  in  the 
satellite  message  which  corresponds  to  a time  mark  in  the  receiver.  Then,  the  counting 
of  the  beat  frequency  begins  at  t + r/c  + e,  where  r is  the  range  at  f . 

b.  In  the  receivers  considered  here,  the  delay  e averages  approximately  1,000  /is. 
By  contrast,  the  delay  in  the  receiver  of  the  received  carrier  is  negligible;  that  is,  the 
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delay  between  the  time  the  signal  is  received  at  the  antenna  and  the  time  it  is  mixed 

with  the  local  reference  frequency  to  produce  the  beat  frequency  is  negligible 

(typically  less  than  10  /js.  This  means  that  the  wave  front  reaching  the  mixer  at  f + At 

+ e is  that  which  was  transmitted  from  the  satellite  at  f + Af  + e — At',  where  At'  is 

the  propagation  delay  of  the  signal  arriving  at  t + At  + e . Again  ignoring  refraction, 

Af'  is  i times  the  range  at  f + Af  + e — At’,  which  is  approximately  r + r(At  + e — 

At'J  . Solving  this  relationship  and  neglecting  second-order  terms  in  £ yields  At'  = Af 

+ -^6,  and  the  time  of  transmission  is  f-i-Af+e  — At'  = f + {1  — — )e  . Thus,  the 
c c 

number  of  cycles  received  at  the  transmitter  between  t,  + Afi  -t-  Ci  and  fj  + Af2  + 

Cj  must  be  equal  to  the  number  transmitted  between  t,  +(1  — rxlc)ei  and  fj  + (1 

- Ic)e2  , or 

f 2 + A f2  +62 

= f,(h~t,  +62-6.  -i^+i^). 

f,  + Af  I + 6 I 

Since  the  exact  delays  arc  unknown  and  cannot  be  determined  from  the  data,  the  last 
two  terms  are  rewritten  in  terms  of  the  mean  delay  e,  the  difference  in  delays  (62  — 
6, ),  and  the  mean  range  rate  r as 


- 6.  6,  ) 


therefore. 


= -«.)+  (r: 

= - Cl ) +^-(>-2  -nW  , 

C C 


t 

I 


,f2  + Af2  +62 


I,  + Af,  + 6, 


f,  + €2  - 6, 


(e?  - e,  ) + i(r2  - r,)e 


] 


c.  The  term  e is  identified  as  the  mean  receiver  delay,  and  may,  on  option,  be 
carried  as  an  unknown  by  the  program.  However,  it  is  indistinguishable  from  a constant 
along-track  or  time  error  in  the  satellite  cphemcris.  and  may  not  be  carried  as  an 
unknown  in  a problem  where  the  satellite  cphemcris  is  also  carried  as  unknown. 
Furthermore,  it  is  also  indistinguishable  from  an  along-track  error  in  the  receiver 
position.  Therefore,  it  cannot  be  solved  for  in  a navigation  solution  which  uses  only  a 
single  pass,  or  even  in  a multipass  solution  using  only  passes  in  a given  direction.  The 
latter  situation  can  arise,  for  instance,  when  a receiver  tracks  only  daytime  passes  of  a 
navigation  satellite  for  a period  of  several  days.  The  restriction  to  only  daytime  passes 
will  cause  the  receiver  to  see  only  the  ascending  (northgoing)  or  only  the  descending 
(southgoing)  branch  of  the  orbit. 
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ti.  When  a reasonably  well-balanced  set  of  northgoing  and  southgoing  passes  is  used, 
the  mean  delay  c is  well  determined.  On  option,  the  program  will  solve  for  a separate 
value  of  e for  each  pass,  for  each  satellite,  or  for  each  receiver.  If  the  receiver  delay  is 
fairly  stable,  it  is  most  reasonable  to  solve  for  a single  value  of  this  parameter  for  each 
receiver  rather  than  for  each  pass.  A separate  solution  would  be  made  for  each  satellite 
if  there  were  a possibility  of  time  biases  in  the  given  satellite  ephemerides.  The  program 
also  allows  a value  for  the  receiver  delay  to  be  input,  in  which  case  the  parameter 
solved  for  is  a correction  to  the  input  value. 


8.  LINEARIZED  OBSERVATION  AND  NORMAL  EQUATIONS,  a.  The  linearized 
observation  equations  arc  obtained  by  expanding  equation  (6)  in  a first-order  Taylor 
series  around  approximate  values  of  the  unknowns  A/'  (frequency  offset).  X 
(receiving  antenna  coordinates),  and  e (mean  receiver  delay).  Thus, 


3AF 


dA 


(7) 


In  the  program,  the  approximate  value  of  AF  at  each  iteration  will  be  cither  32  kHz 
(for  navigation  satellites)  or  16.2  kHz  (for  GEOS  satellites),  although  other  values  may 
be  input  by  the  user.  Initial  approximate  values  of  the  station  coordinates  are  input  by 
the  user  and  improved  at  each  iteration.  The  approximate  value  of  the  delay  c is 
always  zero  on  the  first  iteration. 


b.  In  equation  (7),  d®  is  equation  (6)  evaluated  with  approximate  values  of  the 
unknowns;  6A/',6X  and  6e  are  corrections  to  the  approximate  values.  The  partial 
derivatives  are 


hdc 

dKl- 

= AT  , 

dc/c 

1 

3r, 

dX 

~ X 

“ Jx' 

II 

1 

where 

■ - 

X).  u, 

H 

de 

■ 1 

l). 
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For  the  observation,  let 


^ d -d^ 

c c 


3AF 

be 

H_ 

ax 


Then  the  normal  equation  (for  the  data  from  a single  pass)  can  be  written  as 


A A 12  3 

E" 

aJi  Ai  Bj 

be 

= 

El 

Bi'  C 

- 

F 

where 


Zbjw.b. 


A = Z<ijWa. 

j 

An  = Zaj^jbj 
I 

Ai  = 

B = 

Bj 
C 
E 
E, 

El 


= 

= fm 

T 

= fi  ‘"y"; 

= pw.l^ 

= T.b-W-1- 

J } n 

= TC^W'l- 

Y)  ) 


The  summation  is  taken  over  all  the  observations  during  the  pass.  The  weights,  w-,  are 
equal  to  either  unity  or  to  the  reciprocal  of  the  square  of  the  input  observational 
standard  derivation.  In  addition,  the  quantity 


D = iJ-w-l- 

n n 


is  useful  for  statistical  computations.  In  the  program,  subroutine  NORST  generates  all 
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these  quantities  for  each  pass.  Subroutine  NORSE  is  similar  but  omits  the  calculation 
of  bj  and  the  matrices  .“Ij,  A, 2,  /jj , and  £j  . 

r.  The  three  subroutines,  STASTA,  STASOT,  and  STASOL,  combine  the  matrices 
from  all  the  passes.  They  perform  very  similar  scries  of  matrix  manipulations  designed 
for  a system  of  normal  equations  which,  when  partitioned,  take  the  form 


In  each  case  refers  to  the  pass  unknowns. 

In  STASTA: 

VjIj  = 6A/-  for  the  pass, 

Z = fix  plus  a fie  for  each  satellite  (receiver), 

= A for  the  pass, 

//^  = an  array  made  up  of  .-l,j  and  B for  the  pass, 

= 1:  for  the  pass, 

Q = the  summation  (over  all  passes  of  A2 , and  C,  appropriately 

arranged  by  satellite  (receiver] ), 

R = the  summation  (over  all  passes  of  £’2  and  F,  appropriately  arranged 
by  satellite  [receiver] ). 

In  STASOT: 

= fiAF  and  fie  for  the  pass, 

Z = fix, 

for  the  pass, 

■1,2  A2 
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H,  = ‘ for  the  pass, 

IK  = ” for  the  pass, 

Q = ZC  where  the  summation  is  over  all  passes, 

R = Z/'  where  the  summation  is  over  all  passes. 

In  STASOL: 

= 6AF  for  the  k^^  pass, 

Z = 6X, 

= .-I  for  the  k^^  pass, 

- li  for  the  k*^  pass, 

= E for  the  pass, 

Q = ZC  summation  over  all  passes, 

R = Zf  summation  over  all  passes. 

Also  let 

V = SD  summation  over  all  passes. 

With  a system  of  equations  in  the  form  of  equation  (8),  can  be  eliminated  to  give 

{Q-ZHlG^}H^^)Z  = R-ZHfGj;'Pfc  , 
so  that  the  solution  for  Z is 

Also,  the  sum  of  the  squares  of  the  residuals  is  given  by 

V - ZP[g^‘  Pfc  - (R  - - 2H[Gj*  H^)'*  (R  - 2H[G^‘  P^) 

from  which  their  root  mean  square  (rms)  can  be  calculated.  This,  together  with  the 
inverse  of  the  normal  equations 

(Q-ZW[Gi'Hfc)-' 

corresponding  to  Z,  forms  the  basis  for  calculating  the  error  statistics  for  Z . 
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V'j^  is  calculated  from 


and  the  inverse  of  the  normal  equations  corresponding  to  is 

This,  with  the  rms,  gives  the  basis  for  calculating  the  errors  of  Yj^  . 

9.  WEIGHT  MATRIX  FOR  OBSERVATION  EQUATIONS,  a.  Normally,  each  obser- 
vation is  weighted  inversely  proportional  to  the  variance  of  the  range  difference 
observation.  This  variance  is  defined  in  an  input  array  and  is  assumed  equal  for  all 
observations.  It  is  also  assumed  that  the  range  difference  observations  are  independent, 
having  a random  error  with  zero  mean  and  constant  (but  unknown)  standard  deviation. 
It  IS  possible  within  the  program  to  treat  range  as  the  independent  uncorrclated 
variable  rather  than  range  difference.  This  is  accomplished  by  redefining  the  weight 
matrix. 

b.  Assume  the  range  observation  equation  can  be  written  as 

where  is  the  zero  set  (initial  range)  and  F are  other  parameters.  Further  assume 
that  ranges  (r^.)  are  all  independent  with  equal  uncertainties,  a . 

Then 


and 


ri 
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The  covariance  matrix  of  (rj  — r, , ra  — r, , . . . ^ | is  then 


The  appropriate  weight  matrix  is  then  2'^  and  this  can  be  calculated  using 

p,.i  >.  • _ ^ ^ 

^ j{m  + 1 — i)l(m  + l)o’  for  i>j 

The  introduction  of  this  weight  equation  allows  the  integrated  Doppler  observations  to 
be  treated  as  if  they  were  equivalent  to  range  observations. 

to.  ERROR  ANALY.SIS  CAPABILITY.  An  error  analysis  capability  has  been 
included  in  DOPPLR.  This  capability  permits  the  investigation  of  the  effects  of 
uncstimated  errors  on  estimated  parameters  in  the  Doppler  solutions.  The  effects  of 
biases  and  random  errors  are  calculated  as  perturbations  to  either  observed  or  perfect 
Doppler  counts  and  solutions  performed  using  the  perturbed  data. 

11.  CALPRT  ROUTINE.  The  routine  (CALPRT)  calculates  perturbations  to  the 
observed  Doppler  counts.  Acceptable  error  source  parameters  are  defined  as  follows. 

a Observation  - An  error  in  the  corrected  integrated  Doppler  count.  This  includes 
digitizing,  receiver  phase  shifts,  and  refraction  truncation  errors. 

b.  'I'imin}(— An  error  in  the  time  recovery  within  the  receiver  due  to  varying 
signal-to-noise  ratio.  This  is  calculated  for  each  time  mark. 


.■i 
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c.  offset  I rcijucncy-An  error  in  the  assumed  frequency  difference  between  the 
fundamental  satellite  transmission  frequency  and  the  local  oscillator  frequency. 

d.  Satellite  l-rcjnency- An  error  due  to  an  incorrect  satellite  transmission 
frequency. 

<.  I ropusplierr— An  error  resulting  from  the  limitation  of  the  tropospheric  mixlel. 
This  is  calculated  as  a percentage  error  of  the  modeled  tropospheric  correction. 

/.  Ionosphere- An  error  due  to  neglecting  higher-order  terms  in  using  the 
two-frequency  technique.  This  error  is  modeled  as  a percentage  of  the  measured 
two-frequency  ionospheric  correction. 

g.  hphemens -An  error  due  to  the  uncertainty  in  the  position  of  the  satellite  at  the 
time  of  observation.  The  error  is  subdivided  into  radial,  along-track,  and  crosstrack 
components. 

The  perturbation.  Ad^^  , for  each  error  parameter,  P.  , is  calculated  using 

dd 

• 

I 

where 

d^  - basic  integrated  Doppler  observation  equation  and 
~ the  magnitude  of  the  error  parameter. 

The  calculated  perturbation  may  be  combined  with  the  Doppler  counts  as  a bias,  as  a 
random  perturbation  per  Doppler  count,  or  as  a bias  per  pass  but  random  between 
passes.  The  opticins  available  for  each  parameter  are  shown  in  table  1. 

Table  1.  Parameter  Options 


PARAMETER 

BIA.S 

o/POINT 

0/PASS 

Observation 

X 

X 

Timing 

X 

X 

X 

Offset  Frequency 

X 

X 

Satellite  Frequency 

X 

X 

Troposphere  (%) 

X 

X 

Ionosphere  (%) 

X 

X 

Ephemeris 

Radial 

X 

X 

X 

Qosstrack 

X 

X 

X 

Along  Track 

X 

X 

X 
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SECTION  III.  PROGRAM  OPERATION 


12.  GENERAL  RUN  STRUCTURE.  The  data  deck  necessary  for  a Doppler  solution 
IS  built  from  the  initializing  constants  deck,  the  pass  observation  deck,  and  the  system 
control  cards  as  shown  schematically  in  figure  1.  It  is  possible  to  do  more  than  one 
Doppler  solution  at  a time  by  stacking  decks  similar  to  the  single  solution  deck  as 
shown  in  figure  2,  and  multiple  solutions  using  different  options  with  the  same  set  of 
data  can  be  accomplished  as  shown  in  figure  3. 
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Figure  1.  Single  Doppler  solution  data  deck  structure. 
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Figure  2.  Structure  for  multiple  Doppler  solutions 
using  the  same  options. 


Figure  3.  Structure  for  multiple  Doppler  solutions 
using  different  options. 
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13.  INITlALIZlNd  CONSTANTS  DECK.  a.  The  initializing  constants  which  can  be 
specified  tor  a solution  include  any  of  the  following:  station  coordinates,  spheroid 
parameters,  datum  shifts,  velocity  of  light,  station  receiver  delay,  satellite  frequency 
offset,  satellite  transmitting  wavelength,  title,  and  various  program  control  options.  The 
input  ot  any  of  these  values  is  identified  by  the  appropriate  word  left-justified  on  the 
input  card.  Table  2 lists  the  acceptable  card  formats  to  be  included  m the  deck  of 
constants.  These  inputs  form  the  initializing  constants  deck  and  may  appear  within  this 
deck  in  any  order,  if  more  than  one  card  appears  with  the  same  identification 
left-justified  on  the  card,  the  last  definition  overrides  previous  definitions. 

b.  Input  of  the  perturbations  to  be  included  in  an  error  analysis  run  can  also  be 
included  in  the  initializing  constants  deck.  This  input  requires  two  or  more  cards  in 
sequence;  these  cards,  however,  may  be  in  any  order  with  respect  to  the  other  input  in 
the  initializing  constants  deck.  The  perturbations  are  read  as  a NAMELIST  and, 
therefore,  obey  all  FORTRAN  V rules  pertaining  to  NAMELISTS.  In  addition,  one 
header  card  must  be  included  in  the  NAMELIST  and  all  numerical  input  must  be 
expressed  as  double  precision.  Table  3 gives  acceptable  parameters  in  the  NAMELIST. 
Values  for  any  or  all  of  these  parameters  may  be  specified  in  any  NAMELIST 
definition.  The  random  perturbations  arc  normally  applied  per  observation;  however,  if 
the  value  for  the  random  error  parameter  is  negative,  the  perturbation  is  applied  as  a 
bias  per  pass  but  random  between  passes.  The  card  formats  for  entering  perturbation 
data  are  shown  in  table  4. 

c.  In  computing  a solution  with  real  data,  it  is  only  necessary  to  input  the  OPTION 
and  STATION  cards,  provided  the  station  is  input  on  the  NWL  8D  datum,  if  the 
station  is  input  on  another  datum  specified  by  column  80  on  the  STATION  card,  a 
SPHEROID  card  defining  the  datum  shifts  to  the  geocentric  NWL  8D  datum  is  also 
required.  If  not  input,  the  following  constants  are  preset  to  tlic  indicated  values: 

Title  = solution  for  xxxxxx,  where  xxxxxx  is  the  station  name. 

Velocity  of  Light  = 299  792  500.0  meters  per  second  (m/sec), 

.Satellite  X = 0.749541 21  m (400-1  50  MHz  Mode), 
offset  Frequency  = 32,000  Hz. 

Spheroid  Name  = 3H  NWL, 

AO  = 6 378  145.0  m. 

Flattening  = 0.003352892, 

Dry  Component  of  Surface  Refractivity  = 292.9, 

Wet  Component  of  Surface  Refractivity  = 75.6. 
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Tabic  2.  Deck  of  Constants  Card  Formats 

CARD 


CARD  NAME 

PARAMETER  COLUMN 

FORMAT 

TITLE 

5H  TITLE 

1-5 

A6 

t:ARD 

70-Charactcr  Title 

10-80 

11A6.A4 

STATION 

7H  STATION 

1-7 

A6 

CARD 

Station  Name 

10-15 

A6 

= 0 for  0,X,/i  input  (degrees,  degrees,  meters) 

17 

11 

= 1 for  x,y,:  input  (meters,  meters,  meters) 
Station  Code  (as  defined  on  the  pass  ID) 

18 

A1 

Station  0 or  .v 

20-34 

D15.3 

X or  _y 

35-49 

D15.3 

h or  z (where  h is  the  total  height  of  the  antenna 
above  the  ellipsoid) 

50-64 

D15.3 

Geoid  Separation  (meters)— used  only  in  the  computation  of 
the  tropospheric  correction 

65-79 

D15.3 

Datum  Code- matches  the  first  character  of  the  datum  named 
on  the  SPHEROID  card 

80 

A1 

SPHEROID 

8H  SPHEROID 

1-8 

A6 

CARD 

Datum  Name— first  character  of  this  name  matches  datum  code 
on  the  STATION  card 

10-15 

A6 

Semimajor  Axis  (meters) 

20-34 

D15.3 

Ellipsoid  Flattening 

35-49 

D15.3 

Datum  Shift  (Geocentric  coordinates  of  origin  of  system 
defined  on  this  card) 

X component  (meters) 

50-59 

D10.3 

y component  (meters) 

60-69 

DIO.  3 

z component  (meters) 

70-79 

D10.3 

CONSTANT 

9H  CONSTANTS 

1 9 

A6 

CARD 

Velocity  of  Light  (meters/sec) 

20-34 

D15.5 

Receiver  Delay  (/is) 

35-49 

D1 5.5 

Frequency  Offset  (Hz) 

30-64 

D15.5 

Satellite  Wavelength  (meters) 

65  79 

D15.5 
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CARD  NAME 


Table  2.  Deck  of  Constants  Card  Formats— Continued 

CARD 

PARAMETER  COLUMN  FORMAT 


OPTION 

CARD 


6H  OPTION 


1-7  A6 


OPTION  Punch  Code 

NO.  yes  no 


1 

Number  of  Iterations 

10 

11 

2 

Print  residuals  on  last  iteration  only 

1 

0 

11 

11 

3 

Print  Residuals 

0 

1 

12 

11 

4 

Plot  Residuals 

1 

0 

13 

11 

5 

Observation  Sigma  (meters) 

Preset  to  0.75m  if  not  input 

14-17 

14 

6 

Station  Latitude  Sigma  (meters) 

18-21 

14 

7 

Station  Longitude  Sigma  (meters) 

22-25 

14 

8 

Station  Height  Sigma  (meters) 

26-29 

14 

9 

Plot  station  positions  obtained  from  single 

pass  solutions 

1 

0 

30 

11 

10 

Apply  tropospheric  correction 

0 

1 

31 

11 

11 

Do  single-pass  navigation  solutions 

0 

1 

32 

11 

12 

Data  Source 

33 

11 

= 0 use  real  data 


= 2 use  data  from  previous  run  with 
original  station  input  as  on 
STATION  card 

= 3 use  data  from  previous  run  with 
station  position  from  last  solution 
= 9 Stop 

13  Local  Time  and  Satellite  Track  Plots  46  II 

>1  Print  distribution  of  passes  in  local  time 
= 1 also  plot  subsatellite  points 
= 0 do  neither  of  above 


14  Perform  solution  without  time  bias 

adjustment  0 1 47  11 

15  Delete  samples  with  elevation  angles  less  than 

this  value  (degrees)  48-49  12 
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Tabic  2.  Deck  of  Constants  Card  Formats -Continued 

CARD 

CARD  NAME  PARAMETER  COLUMN  FORMAT 


OPTION  Punch  Code 


NO. 

yes 

no 

16 

Local  Time  Deletion  of  Passes 

= 0 no  deletion  on  local  time 
= XX  and 

50-51 

12 

= yy  keep  those  with  xx<local  time<^y 

52-53 

12 

17 

rms  multiple  for  residual  editing 

54 

11 

18 

North— South  Pass  Deletion 

55 

11 

= 0 use  all  passes 
= 1 use  only  north-going  passes 
= 2 use  only  south-going  passes 

19 

Ecjuipment  Delay  Adjustment 

56 

11 

= 1 adjust  e varying  pass  to  pass 
= 2 satellite  to  satellite 

= 3 receiver  to  receiver 

20 

Apply  ionospheric  correction 

0 

1 

57 

11 

21 

Fix  frequency  offset  for  single-pass 

solution  plots  to  values  from  previous  • 
solution 

2 

0 

58 

11 

22 

Deletion  multiple  for  single-pass  solution 

plots 

59 

n 

23 

I3ata  Type 

60 

11 

= 0 nr  Doppler  data  {P,  T,  R) 

= 1 Quick-check  data 
= 2 Backpack  data 
= 3 Geocciver  data 
= 4 ITT  Doppler  data  (P,  7'^) 

24 

Timing  Mode 

61 

11 

= 0 satellite  clock 
= 1 local  clock 

25 

Generate  Perfect  Data 

9 

0 

62 

11 

26 

Weight  Matrix 

63 

11 

= 0 range  difference  formulation 
= 1 range  formulation 
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Table  3.  NAMELIST  Parameters 
ERROR  SOURCE  ACCfin'AlJl.F.  PARAMETER  NAMES 

HI  AS  RANIX)M  ERROR 


Observation 

OBIAS 

OBSIG 

Timing 

TIMEB 

TIMSIG 

Offset  Frequency 

FBIAS 

fsk; 

.Satellite  Frequency 

— 

tfsk; 

Tropospheric 

TROPOZ 

— 

Ionospheric 

Ephemeris 

lONOZ 

Radial 

RADIAL 

SIGRAD 

Crosstrack 

CROSS 

SIGCRS 

Along  Track 

ALONG 

SIGALG 

Table  4. 

Perturbation  Data  Card  Format 

CARD 

CARD  NAME 

PARAMETER 

COLUMN  FORMAT 

CARD  1 

5H  SPERT 

1-5 

A5 

CARD  2 

.SH  SPERT 

1-5 

A5 

CARDN 

4H  SEND 

1 4 

A4 

14.  OB.SERVATION  DATA  DECK.  The  observation  data  deck  for  each  pass  consists 
of  a pass  identdication  (ID)  card,  ephcnieris  cards,  and  Doppler  cards  (figure  4). 


ICAVATNMOATA  MCC 


Figure  4.  Observation  data  deck  structure. 

a.  I’lisa  liicHtijicatiou  Card.  The  pass  ID  card  is  the  first  card  for  each  pass  and  must 
follow  the  format  showfn  in  table  5.  The  quantities  in  columns  1—26  are  required 
inputs.  If  the  wet  and  dry  bulb  temperatures  and  pressure  are  not  specified,  they  are 
set  to  standard  values.  On  the  pass  ID  card,  if  either  of  the  values  for  frequenev  offset 
or  satellite  wavelength  is  zero,  preset  values  or  values  from  the  constants  deck  are  used. 
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Table  5.  Pass  Identification  Card  Format 


CARP  NAME 


CARP 

PARAMETER  COLUMN  FORMAT 


PASS  ID 


4H  PASS 
Pass  Number 
Satellite  Identification 
Station  Identification 
Montii  of  Pass 
Day  of  Pass  (GMT) 

Year  of  Pass 
Hour  of  Pass  (GMT) 

Minute  of  Pass  (GMT) 

Mode  (Gcoceiver  input  only) 
I>y  Bulb  Temperature  (°F) 
Wet  Bulb  Temperature  (°F) 
Pressure  (mm/Hg) 

I>v  Bulb  Temperature  (°C) 
Relative  Humidity  (%) 

Pressure  (mb) 

Frequency  Oflset  (Hz)  tor  pass 
Satellite  Wavelength  (m)  for  pa 


1 -4 

A4 

5-10 

A 6 

11 

A1 

12 

A1 

14-15 

12 

16-17 

12 

18-19 

12 

23-24 

12 

25-26 

12 

28 

12 

For 

31-40 

D10.2 

observation  types 

41-50 

D10.2 

0,  1,  and  2 

51-60 

D10.2 

For 

31-40 

DIO. 2 

observation  types 

41-50 

D10.2 

3 and  4 

51-60 

D10.2 

61-70 

D10.2 

s 

71-80 

D10.2 

h.  C.lock  l-pocli  Cjxljte  Ciird.  The  epoch  update  card  permits  a local  clock  epoch 
error,  drift  rate,  and  scale  factor  to  be  introduced  in  the  preprocessor.  This  card,  if 
used,  follows  the  pass  ID  card  and  all  subsequent  observational  times  arc  corrected 
using  this  information  until  another  epoch  update  is  encountered.  This  input  can  only 
be  used  witli  HT  and  ( ieoceiver  observations.  The  card  format  is  shown  in  table  6. 

Table  6.  Fpoeh  Update  (.^rd  Format 

. CARD 

CARD  NAME PARAMETER  COLUMN  FORMAT 


ep(k:h  update 

6HTIME 

1 -4 

A4 

CARD 

DAYS  since  Jan  0 

10-12 

13 

Hours 

15-16 

12 

Minutes 

18-19 

12 

Epoch  Delay  (/is) 

25-39 

D15.10 

I'lrift  Rate  (^s/day) 

45-59 

D15.10 

Scale  I'actor 

65-79 

D15.10 

.R) 


c.  t.i>ln  iiivris  i'.ard  Input.  The  satellite  eplieineris  for  each  pass  can  be  input  to  the 
pr.  )cessor  111  any  of  the  three  following  formats: 

• As  discrete  satellite  ephemeral  points  at  a fixed  time  interval. 

• As  coefficients  lor  Chebychev  polynomials  which  describe  the  fixed-Earth  x,  y, 
and  .r  variation  of  the  satellite  with  time.  7’hesc-  arc  produced  by  the  separate  program 
COEYTP  from  an  ephemeris  on  magnetic  tape. 

• As  .satellite  mes.sage  words  received  from  the  navigation  satellite  (predicted 
ephemeris). 

(1)  Ephemeris  Point  Input,  if  the  ephemeris  is  input  as  discrete  .satellite  points, 
the  structure  of  the  ephemeris  deck  follows  the  format  shown  in  table  7.  Following 
card  1.  the  satellite  points  are  input  one  per  card,  as  indicated.  The  card  for  the  last 
ephemeral  point  must  have  a “1”  in  the  first  column  to  signify  the  end  of  tlie 
ephemeris  input.  The  tunc  for  each  ephemeris  point  must  increase  by  N2  minutes 


Table  7.  Ephemeris  C^ard  Format  (Discrete  Satellite  Points) 


CARD  NAME 

HARAMETER 

CARD 

COl.UMN 

FORMAT 

CARD  1 

6H  POINTS 
4H  PASS 

I r> 

2f  24 

A6 

tlrbit  Number 

27  32 

A6 

EPHEMERIS 

Satellite  Identification 

33 

A1 

POINT  CARD 

Termination  Indicator 
= 0 normally 
= I last  ephemeral  point 

1 

11 

Month  of  Ephemeral  Point  (CMT) 

9-10 

12 

Day  of  Ephemeral  Point  (CMT) 

11  12 

12 

Year  of  Ephemeral  Point  (CJMT) 

13-14 

12 

Hour  of  Ephemeral  Point  (CMT) 

16  17 

12 

Minute  of  Ephemeral  Point  (CMT) 

Pixed  Earth  Ephemeral  Position  and  Vel 

18-19 

ocity 

12 

A'  (meters) 

2 1 - 30 

FI0.3 

y (meters) 

31-40 

FI  0.3 

' (meters) 

41-50 

F10.3 

.V  (mcters/sec) 

5 1 60 

F10.3 

y (mcters/sec) 

61-70 

F10.3 

; (mcters/sec) 

71  80 

FI  0.3 

31 


over  the  preceding  point.  N2  is  a parameter  pre.set  to  2 minutes  {preprocessors  PR(X^ 
or  OBSIN)  and  corresponds  to  the  integration  interval  over  which  the  Doppler  counts 
arc  measured.  A missing  ephemeris  point  will  cause  the  entire  pass  to  be  deleted.  The 
maximum  number  of  ephemeris  points  which  can  be  entered  for  a particular  pass  is 
controlled  by  parameter  N1  which  is  preset  to  1 1.  An  example  of  the  input  card  deck 
for  the  ephemeris  points  is  shown  in  figure  5.  These  points  are  for  pass  13504  of 
satellite  Y on  23  June  1970.  The  first  ephemeral  point  is  at  05  hours  16  minutes  and 
the  last  at  05  hours  36  minutes. 


(2)  Ephemeris  Coefficient  Input.  The  segment  of  the  ephemeris  over  which 
observations  have  been  made  can  be  fit  with  Chebychev  polynomials  of  the  form 


n 

n + rn 

x!(N  — w)! 

tn 

rn 

(x  - m)!N! 

where  the  observation  interval  has  been  divided  into  N equal  segments  so  that  x 
takes  the  values  (0,  1,  2,  ...,  N).  The  coefficients  can  then  be  input  to  the 

processor,  thus  reducing  the  input  needed  to  specify  an  ephemeris  arc.  The  formats  for 
the  ephemeris  deck  cards  with  Chebychev  coefficients  are  shown  in  table  8. 

An  example  of  this  ephemeris  input  for  pass  13504  is  shown  in  figure  6.  From  these 
coefficients  the  actual  ephemeral  points  can  be  evaluated  by  the  following  formula; 
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t = 0 
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e 6.  Listing  for  cphcmeris  coefficients  input. 


The  coefficients  shown  in  figure  6 were  computed  by  the  ephcmeris  preprocessor  from 
the  1 1 points  indicated  for  this  pass  in  figure  5.  These  points  are  obtained  from  the 
ephemeris  tape  by  the  ephemeris  preprocessor,  which  searches  the  ephemeris  tape  and 
either  extracts  the  discrete  points  or  computes  the  coefficients  and  outputs  the 
ephemeral  segment  on  punched  cards  with  the  formats  shown. 


Table  8.  Ephemeris  Card  Format  (Chebychev  Coefficients) 

CARD 

CARD  NAME  PARAMETER  COLUMN  FORMAT 


CARD  1 

12H  COEFFICIENTS 

1-12 

A6 

4H  PASS 

18-21 

— 

Orbit  Number 

23-28 

A6 

Satellite  Identification 

29 

A1 

6H  POINTS 

33-38 



Number  of  points  used  in  Chebychev 
polynomial  fit  for  each  component 


/^./^.aiid^ 

40-41 

12 

6H  DEGREE 

47-52 

— 

Degree  of  Polynomial  Fit 

55 

— 

4HTIME 

61-64 

— 

Time  interval  in  minutes  between 

satellite  points  used  in  polynomial 
fit 

66-67 

12 

CARDS  2-  4 IHX  or  IHY  or  1 HZ  corresponding  to 

input  of  the  polynomial  coefficient 

for  the  X component,  y component, 
or  ; component  of  the  ephemeral  arc 

1 

A1 

Hour  of  first  point  used  in  polynomial 
fit  wliich  corresponds  to  (.v  = 0) 

(GMT) 

7-8 

12 

Minute  of  first  point  (GMT) 

9-10 

12 

'3x0 

11-20 

D10.4 

■3x1 

21-30 

D10.4 

•3x  j 

31-40 

D10.4 

3x3 

41-50 

D1 0.4 

3x« 

51-60 

D10.4 

■3x5 

61-70 

D10.4 

-dx* 

71-80 

D10.4 
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(3)  Satellite  Message  Eplienieris  Input. The  iliird  form  of  cplicmeral  input  is  the 
satellite  message  ephemeris  as  received  trom  NNS.  The  satellite  message  consists  of 
both  fixed  parameters  and  variable  parameters;  a complete  description  of  these 
parameters  can  be  found  in  Gutheim.  Table  9 gives  the  card  format  for  the  satellite 
message  ephemeris  input.  The  parameters  entered  on  cards  2 through  5 are  exactly  as 
received;  the  fixed  parameters  are  decoded  in  the  Doppler  processor  into  orbital 
elements  and  the  variable  parameters  into  2-minute  corrections  to  the  orbital 
parameters.  The  number  of  variable  parameters  entered  on  cards  2 and  3 is  determined 
bv  the  timespan  of  tiie  observations  for  the  particular  pass.  The  first  variable  parameter 
must  precede  the  time  of  the  first  observation  and  the  last  variable  parameter  must  be 
later  than  the  last  observation  time.  If  less  than  nine  variable  parameters  are  input,  a 
blank  card  must  be  used  for  card  3. 

t/.  Doppler  Input  Deck.  Three  different  formats  can  be  used  for  the  input  of  the 
integrated  Doppler  observations.  The  type  of  Doppler  input  is  specified  in  the  option 
card  (column  60)  and  must  be  one  of  the  following: 

• Geoceiver  • ITr  5500  • Backpack. 

Within  a solution  deck,  the  Doppler  counts  must  be  of  a single  type;  however,  in  multiple 
processing  runs,  the  Doppler  data  type  may  differ  between  successive  solution  decks. 

(1)  Geoceiver  Doppler  Input.  A Geoceiver  observation  consists  of  a time 
indicating  the  start  of  the  counting  interval,  the  duration  time  of  the  count,  the 
Doppler  count,  and  the  refraction  count.  Two  observations  may  be  put  on  a single  card 
in  the  format  shown  in  table  10.  The  Doppler  deck  consists  of  observation  cards  only. 
The  time  sec]uence  of  the  observations  on  the  cards  or  within  the  deck  is  not  important 
in  the  program  executum  but  mav,  however,  result  in  a difference  in  internal  storage 
requirements.  Internal  program  parameters  limit  the  number  of  ephemeris  points  to  35 
per  pass  and  these  points  arc  controlled  by  both  the  number  of  Doppler  counts  input 
and  their  time  sequence. 

(2)  ITT  5500  Doppler  Input.  The  input  of  ITT  5500  integrated  Doppler 
obsc-rvations  is  identical  to  the  GeiKciver  input  except  that  the  input  of  seconds  for  the 
starting  time  and  duration  time  is  rcpl.aced  by  line  counts.  One  line  count  is  equal  to 
(23-4/6103)  ' 120  seconds.  This  simplifies  the  input  of  the  ITT  Doppler  counts  since 
this  unit  outputs  a Doppler  count  and  refraction  count  every  4.6010159  seconds  as  it 
is  accumulating  the  full  2-mmute  count.  A special  symbol,  “P,"  which  is  used  in  the 
column  following  the  duration  line  count,  indicates  that  an  interval  of  0.3735868 
seconds  has  been  included  in  the  observation  duration  (i.e.,  26  lines  -t-  0.3735868 
seconds  = 1 20.0  seconds).  The  format  for  input  of  ITT  observations  appears  in  table  1 1 . 

* CJUTHEIM.  (.  C.  “Program  Roquircmcno  for  Two  Minutr  Integrated  Doppler  Satellite  Navigation  Solution.’ 
\pptUd  Phytui  LMht>ratory  I'rchnifal  Mrmoramium  TG  819  1.  Silver  Spring,  Md.:  Johns  Hopkins  Universitv  . 
April  1967 
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CARD  NAME 
CARD  1 

CARD  2 
CARD  3 
CARD  4 
CARD  5 


Table  9.  Satellite  Message  Card  Format 


PARAMETER 

CARD 

COLUMN 

FORMAT 

9H  SATELLITE 

1-9 

A6 

Orbit  Number 

21-26 

— 

Satellite  Identification 

27 

— 

Variable  Parameter 

1 

1-10 

no 

2 

11-20 

110 

3 

21-30 

110 

4 

31-40 

no 

5 

41-50 

no 

6 

51-60 

no 

7 

61-70 

no 

8 

71-80 

no 

Variable  Parameter 

9 

1-10 

no 

10 

11-20 

no 

1 1 

21-30 

no 

12 

31-40 

no 

13 

41-50 

no 

14 

51-60 

no 

15 

61-70 

no 

16 

71-80 

no 

Fixed  Parameter  Inputs 

Time  of  Perigee 

1 -10 

no 

Mean  Motion  N 

11-20 

no 

Argument  of  Perigee  at 

21-30 

no 

Precession  Rate  of  Perigee 

31-40 

no 

Eccentricity 

41  50 

no 

Mean  Semimajor  Axis 

51-60 

no 

Right  Ascension  of  Ascending  Node  at 

61-70 

no 

Precession  Rate  of  Node 

71-80 

no 

Fixed  Parameter  Inputs 

Cosine  of  inclination 

1 10 

no 

Ixjngitude  of  Greenwich  at 

11-20 

no 

Time  of  Last  Satellite  Injection 

21-30 

no 

Sine  of  Inclination 

41-50 

no 

Frequency  Offset  in  ppm 

51-60 

no 
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CARD  NAMl 

T.ible  10.  Doppler  ( i.ird  Fisrmat  (Ceoceiver) 

lARI) 

PAKAMFTFR  COLUMN 

FORMAT 

OB.SERVA  riON 

Orbii  Niiniber 

1 6 

A 6 

CARD 

Satellite  Identification 

7 

A1 

Station  Identification 

Starting  Time  ((JMT)  for  ( )bservation  1 

8 

A1 

Hour 

10  II 

12 

Minute^ 

12  13 

12 

Seconds 

Counting  Duration 

14  19 

F6.4 

Minutes 

20  21 

12 

•Seconds 

22  30 

F9.6 

Doppler  Count  (decimal) 

3 1 38 

18 

Refraction  Count  (decimal  + 2000) 
.Starting  Time  (CMT)  for  ( )bservation  2 

39  43 

15 

Hour 

45  4f) 

12 

Minutes 

47-48 

12 

.Seconds 

Counting  1 )uration 

49  54 

F6.4 

Minutes 

55  56 

12 

.Seconds 

57  65 

F9.6 

Diqipler  Count  (decim.il) 

66  73 

18 

Refraction  Count  (decimal  + 2000) 

74  78 

15 

(3)  li.iikpack  Dop|ilfr  Input.  I’lie  Poppli-r  prcpmtcssor  for  the  Backpack  data 
(PRO(')  retjuircs  only  an  e.itimate  of  the  starting  time  ot  the  first  Doppler  count  in 
order  to  m.itch  tlie  IXippler  counts  to  the  correct  satellite  2-minute  interval.  All  the 
Doppler  ctiunts  must  he  of  2-mmute  duration  and  must  be  corrected  for  ionospheric 
refraction.  I'able  12  gives  the  input  format.  The  program  parameters  are  set  to  allow 
two  Doppler  cards  at  most  or  eight  2 minute  Doppler  observations;  however,  for  four 
or  less  observations  onl\  one  card  need  be  input. 

15.  OL'TPl  r DKS('KII’'MON.  I’he  output  from  the  DOI’l'l.R  program  is  controlled 
by  the  options  specifieil  on  the  OPTION  caul  in  the  initiali/ing  lonstants  deck.  On 
initial  entr\  into  the  program,  all  input  data  are  read  and  a data  file  is  created  by  the 
preprocessor.  During  this  process,  the  output  shown  in  figure  7 is  generated.  The 
initializing  const.ints  for  the  particular  run  are  output  along  with  a pass-by  pass 
summary  ot  the  data  and  associated  corrections.  Following  the  generation  of  the  data 
file,  the  specific  data  processing  options  on  the  OIH'ION  card  are  executed  m sequence. 
A description  cjf  each  of  the  data  piocessing  options  and  the  associated  output 
generated  by  the  options  are  as  follows. 
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CAKl)  NAME 


Table  1 1.  [Soppier  Card  Format  (ITT) 

CARD 

PARAMETER  COLUMN  FORMAT 


OBSERVATION 

CARD 


Orbit  Number 

1-6 

A6 

.Satellite  identification 

7 

A1 

Station  Identification 

8 

A1 

Starting  Time  (GMT)  for  Observation  1 

Hour 

10-11 

12 

Minutes 

12-13 

12 

Line  Counts 

14-15 

12 

Qjunting  Duration 

Minutes 

20-21 

12 

Line  Counts 

22-23 

12 

Special  Interval  Symbol  “P” 

24 

A1 

Doppler  Count 

31-38 

18 

Refraction  Count 

39-43 

15 

Starting  Time  (GMT)  for  Observation  2 

Hour 

45-46 

12 

Minutes 

47-48 

12 

lane  Counts 

49-50 

12 

Counting  Duration 

Minutes 

55-56 

12 

Line  Counts 

57-58 

12 

Special  Interval  Symbol  “P” 

59 

A1 

Doppler  Count 

66-73 

18 

Refraction  Count 

74-78 

15 

cl.  <)l’'ll().\  9 /’/of  Suition  /’ovifionx  Obtained  I roin  Sini^lc-I'ass  Solutions.  Three 
single  pass  navigation  solutions  arc  performed  constraining  tlic  station  height, 
longitude,  and  latitude,  respectively,  for  each  pass.  The  resulting  change  in  station 
position  referenced  to  the  input  coordinate  is  output  as  shown  in  figure  8.  Plots 
(figures  9.  10,  and  II)  are  also  generated  for  each  of  the  three  different  constraint 
conditions  showing  the  distribution  of  the  navigation  solutions  for  each  pass.  A 
deletion  rms  multiplier  (OPTION  22)  can  be  input  so  that  extreme  navigation  values 
will  not  be  plotted.  In  addition,  if  a combined  solution  was  performed  prior  to  the 
navigation,  OPTION  21  will  also  constrain  the  value  of  the  frequency  offset  for  each 
pass  used  in  the  navigations  to  that  ol  the  multipass  .solution. 

h.  ()l’'ri()\  It  Do  Si nf^lc-IKiss  iWieifitition  Solutions.  Single-pass  navigation  solu- 
tions are  computed  in  the  plane  defined  by  the  satellite  positions  corresponding  to  the 
first  and  last  observation  and  the  input  station  position.  The  two  components  of  the 
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Tabic  12.  Doppler  Card  Format  (Backpack) 

CARIi 


CARD  NAMK 

PARAMF.TF.R 

COl.UMN 

FORMAT 

OBSERVATION 

Orbit  Number 

1-6 

A6 

CARD 

Satellite  Identification 

7 

A1 

Station  Identification 

8 

A1 

Station  Name 

9 14 

A6 

Doppler  Type  Code 
= 1 Data  = Obs 
= 2 Data  = Obs 
Observation 

18 

11 

1 

20-30 

01 1 

2 

, For  data  type  2 

35-45 

on 

50-60 

01 1 

4 

Observation 

65-75 

on 

1 

24-30 

17 

2 , 

39-45 

17 

^ For  data  type  1 

54-60 

17 

4 

69-75 

17 

station  position  in  this  plane  and  a frequency  offset  arc  estimated  for  each  navigation 
solution.  A total  of  three  iterations  arc  performed  with  residual  editing  based  on  the 
rms  times  a multiple  which  has  been  specified  by  OPTION  17.  The  single-pass 
navigation  output  is  shown  in  figure  12. 

c.  0/''/70,\’  14  Perform  Solntio>i  U'ilhoiil  Time  ULix  Adjuxtmetit  A multipass 
solution  is  performed  in  which  the  estimated  parameters  are  the  three  components  of 
the  station  position  and  a frequency  offset  for  each  pass.  The  output  for  this  option  is 
identical  to  the  output  generated  by  Of'TION  19,  except  item  4 does  not  appear. 

(/.  ()l’'no\  19-Iujuipmcnt  Delay  Ailiustments.  A multipass  solution  is  performed 
m which  the  estimated  parameters  are  the  three  components  of  the  station  position,  a 
frequency  offset  for  each  pass,  and  one  or  more  equipment  delay  parameters 
(paragraph  7).  Provisions  are  made  in  the  specification  of  this  option  which  enable 
equipment  delay  parameters  to  be  defined  for  either 

• Fateh  pass, 

• Each  satellite,  or 

• Each  receiver. 
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LA.impIcs  lor  tlu'  output  of  c;u  li  ot  tlu'st-  arc  sliowii  in  figures  13  aiul  14.  For  tlic 
''.ticllitc  or  receiver,  an  cc]uipmeiu  ciel.iv  p.iramcter  is  estimated  lor  eaeli  different 
satellite  or  station  code  (pass  identification  card)  found  in  tlie  observation  data  deck 
used  in  tlie  nuiltip.iss  solution. 

f.  Ol'IKiS  4 I'lot  licsiiliiiils  The  observation  residuals  from  a multipass  solution 
e.in  be  [ilotted  usinr*  this  option.  The  output  generated  for  a single  pass  is  shown  in 
I'lgurc  1 3. 

/.  orrioV  /)  I’iisf  Distiilmtion  I’lol!:.  The  distribution  of  passes  with  respect  to 
local  st.ition  time  and  with  respect  to  the  tracking  station  (tigure  1(>)  can  be  generated 
with  this  option. 

1(>.  Ol  TIH  r l-KROR  MFSSA(.LS.  Error  messages  produced  In  the  preprocessors 
OBSIN  imbor  I’lUK'  are  as  follows: 

.1  ALL  PASSES  DEI.ETEI) RUN  TP.KMINATED  No  passes  were  output  on 

drum  bv  the  preprocessor  PRCKl. 

h DOPPLER  AND  EPHEMERIS  ARE  NOT  CONSI.STENT  PASS  DELETED 
Epheiiiens  points  input  to  the  preprocessor  PRCK)  are  not  increasing  at  2 nimute 
intervals. 

, 1X)PPLERS  ARE  MISSINCI-Doppler  counts  aic  missing  for  the  cpiick  check  data 
type. 

./  EPHEMERIS  COEFFICIENTS  MISSINC  E(3R  THIS  PASS-  The  ephemeris  data 
deck  is  mwsmg  tor  this  pass. 

I.  F.KROK  IN  EPHEMERIS  INPI.IT  An  er'or  w.is  detected  in  c.irds  2 to  5 of  the 
ephemeris  c oc'ffii  lent  input.  The  coefficients  were  either  missing  or  mispunched. 

/ l:RKOK  IN  INPUT  PASS  DELETED- An  error  w.is  detected  wliile 

preproc  esMiig  the  ephemeris  input  deck  or  Doppler  input  deck.  Tlu  pass  was  deleted 
troni  liiither  processing. 

ti  INPUT  DATA  SKIPPED  Input  dat.i  cards  were  skipped  in  search  tor  next  pass 
observ.itioii  deck. 

/(  INPU'T  ERROR  An  error  was  detectol  while  preprocessing  the  initiali/iiig 
const.ints  dec  k.  The  card  giving  the  error  is  printed. 
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I l.BSS  rHAN  XX  DATA  POINFS  AITFn<  SEQUENCMNC  PASS 

Of'.I.ElEI)  rills  error  octuretl  iii  tlie  preprocessor  PK(K'  wliile  proeessins;  cither 
lhiekp.iek  or  ijmek  clu\k  il.it.i.  I here  were  less  th;m  XX  data  points  after  seijiienciiig 
the  observations  in  iiiireasin^  order  prior  to  niatehing  them  to  the  correct  ephemeris 
intervals. 

) MIMPER  ()l  SAMPI.es  EXCEEDS  MAX  ARRAY  LENCTU  The  number  of 
ephemeris  points  evaluateii  for  this  jiass  esceederl  ,35,  which  is  a preset  program 
parameter  lor  the  mavinuim  satellite  points  per  pass.  Additional  Doppler  observations 
for  this  pass  are  not  usi'd. 

k.  PASS  XXXXXXXX  fS  FOR  ANOTHER  STA I ION  The  station  identification 
on  the  PA.SS  identification  card  does  not  agree  with  the  station  identification  on  the 
STATION  input  canf. 

/ PREDICTED  AND  MEASURED  DOPPI.ERS  COULD  NOT  BE  MATCHED-The 
preprocessor  PR(X  could  not  match  the  observed  Doppler  counts  to  the  predicted 
Dojipler  counts. 

m.  RUN  MISSIN(;  EITHER  S I ATION  OR  SPHEROID  INPUT-Tliis  message  results 
if  the  S PA  i ION  tir  ORPION  cards  aie  missing  from  the  initializing  constants  deck  or  if 
the  STATION  card  specifies  a datum  type  other  than  NWL  and  the  corresponding 
SPHEROID  card  is  missing. 
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Figure  9.  Single-pass  navigation  solution  plot-height  constrained. 
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residuals  repeated  for  each  pass 
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RESIDUALS  REPEATED  FOR  EACH  PASS 


r 


© 


% 

< 

a. 

I 

O 

< 

UJ 

oc 

O 


Q 

LU 

H 

< 

Ui 

a. 

lij 

X 

CO 

< 

3 

Q 

CO 


0 


_|  M Cj  C>  o 


— to  o 

o 

— M C3 


o-  o lO 
< o « 
J o o « 


o • — 
o < 

3 0 — 


»■  — 
in  ♦ 

• n« 


X O X ► 

■ a V « 

1 — j 

« h.’  lu  r>  — C* 

^ X » • C O 

e —no 


— — c- 


C — IT 
O n*  — 


in  n 

> • • • • 
• w • • 


• ♦ o — 
soon 
»-  « o c> 


U)  o ♦ ^ 
19  O O • 
O « — 


ft  n> 

o « 
or  o 
>■  ft  • 


a 7 a 
•*  O 0 

3 — • 

n *-  « 
u u n 


ft  0 « 0 — 


K « -e  o r>  0 r« 

ft  0 ft  — a — rn 
ft  V ft  m u>  m in 
oo  o o o o o 
nm  »n  r»  pn  r»  « 
lftP4  fft  CM  Pft  P*«  fft 
ft  ft  ft  ft  ft  ft  ft 


ft  0 ft  — O — n« 

ft  « ft  in  in  in  in 
oo  0 3 0 0 0 
r^f<»  M M m n 


Oo  o o o o o 


k.7  OO  O O O O O 

cnnoccrt 
c c o C'  n o c*  c- 

•—  On  o o o o n 


O—  ■ftnftC'  — 
ft  ft  ft  in  in  ft  « 


I PM  n«  pft 


• ft 

111  0 
r • 


► 7 

C ft 

»>  1 


ft  *•  n 

o • o 

n*  ft  O 

C n#  Pm 


i 

* 

(O 

CM 

? 


O 

h- 

a. 

O 


£ 


II 

5 « 


ll 

> o 
c •= 

i c 
= ° 

& f 

£ a 

S 3 
= E 
s -g 


° I 

CM  g 

^ a. 

o ; 


J 8 

o “ 
”■  c 
2 2 


II 


»_  o 
— ^ 


s r 


; 8 
V — 
X 9 

c a 


tr  .2 
« “ 
m 


• t _; 

■2  I § 

ft  ^ o 

I I i 

& 'a  Al  £ 

1 i ? I 

E c g e 

S 2 9 t 

I P 5 ? 

S o 2 ■§ 


5 5 ^ 


3 2 e 5 


•2  £ r. 


a>  — 


i «j  c 

E O'  ~ 


2 ft 


o 

!>• 

"o 


o 

k. 

T3 


>> 


3 

cr 

o 

-r 

■5 

c 

0 

«-> 

o 


E ;5 

£ a = 
.1?  £ ^ 


s 


z 

0 

h- 

c 

*0 

0 

0 

5 

? 

t) 

JZ 

1 

k 

g 

ft 

0 

1 

0 

>• 

igure 

cu 

0 

c 

g 

*o 

1 

ft 

u 

& 

u. 

4) 

ft 

■s 

c 

0 

u 

C 

! 

3 

? 

ft 

c 

i 

c 

a 

ft 

ft 

JZ 

1 

u 

c 

0 

s 

1 

c 

S 

c 

0 

1 

0 

0 © 0G00O 


51 


52 


Pass  code  for  Mio:iite  distribution  plot 

Indication  that  pass  22861 YF  was  a northsoina  to  vaulhgoing  pass  which  occurred  between  0200  and  0300  tours  local  time. 
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sr;c  NON  IV.  i;pmi  mi  ris  PKn»K(K:r.ssoR 

17.  (iENERAL  l)ESCRIPTK)N.  a The  eplitmcris  prepmccssor  is  tlcsignccl  to  process 
the  epliemeris  tapes  ereateii  b\  the  IVIeiise  Mapping  Agency  Topograpli ic  Center. 
(These  ephementles  were  previously  produced  by  tiie  U.S.  Naval  Surface  Weapons 
Center  | NSWC]  and  follow  the  format  .is  described  by  Sims  * ).  Each  tape  contains  tlic 
position  and  velocity  components  for  one  NNS  at  I minute  UTC  time  intervals.  A 
single  tape  may  contain  this  ephemeral  mform.ition  for  a period  of  up  to  2 weeks.  To 
simplify  the  input  to  the  Doppler  processing  program.  COPYTP  searches  the  ephemeris 
tape  and  extracts  segments  of  the  epheiiu  ris  coineident  with  the  Doppler  observations. 
The  segments  to  be  extracted  are  specified  b\  input  data  cards  to  fXlPYTP. 

h.  Chebychev  poKnomials  .iri-  fit  to  each  segment  extracted  from  the  ephemeris 
tape  thus  allowing  interpolation  by  the  |rroeessor.  The  rms  of  the  residuals  of  the 
polynomial  fit  is  calculated  and.  if  it  is  within  a predetermined  tolerance,  the 
polynomial  coefficients  are  output  from  (.'OPYTP.  If  the  poivnomials  do  not  fit  the 
ephemeral  arc,  the  actual  eph  •mens  points  are  output.  Either  type  of  output  is 
acceptable  as  input  to  the  Doppler  processor. 

18.  EPIIEMI'^RIS  I'AI’I.  I'ORMAT.  The  ephemeris  tape  as  supplied  by  NSWC  has  the 
position  and  velocity  components  for  one  NN.S  at  1 -minute  time  intervals.  The 
ephemeris  is  m an  Earth  fixed  geocentric  geodetic  coordinate  system.  One  record  of 

i the  binary  coded  decini.il  (BCD)  card  ini.ige  tape  is  formatted  as  shown  in  table  13. 


NAMI-. 

NWl.  TAI’E 


I 


'I'able  13.  BCD  Card  linage  Tape  Format 

CAKI) 


l•AKAM^^;l< 

COLUMN 

LOKMAT 

Last  two  digits  of  year 

1 

-) 

12 

Epoch  (days  since  J.in  0) 

4 

(> 

13 

Hours  (CMT) 

H 

') 

12 

Minutes  (CMT) 

1 1 

12 

12 

.V  eomponeiit  of  position  (km) 

1 4 

22 

FI  0.3 

y com|ronent  of  position  (km) 

2 1 

42 

F10.3 

e eoin|ronent  of  position  (km) 

43 

42 

F10.3 

.V  component  of  velocity  (km, -see) 

4 4 

52 

FI0.4 

y component  of  velocitv  (km/see) 

S3 

fi2 

FI  0.4 

c component  of  veloc  itv  (km/see) 

f)4 

72 

FI  0.4 

.Satellite  Number 

77 

80 

14 

* SIMS.  T.  "NWL  IVcti^inti  Ephemeris."  Saval  lv'rup*»ns  l.ahitriitorv  Trchnual  Rt  junt 
N4V4I  Surfaic  Weapons  Onicr.  September  1972. 


I'K  >842  ILihluTen.  Va.: 
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\ 19.  DATA  DECK  STRUCTURE,  a.  liplwmvris  Preprocessor  Run  Deck.  The  run 

deck  structure  for  the  ephemeris  preprocessor  COPYTP  is  shown  in  figure  17.  PROG  is 
the  program  tape  and  DAP26  is  tlie  N.SWC.'  epliemeris  tape. 


r»*OTCO?VTf 

I •CO^IM.A  WK>G  CO^VTl^ 
j IM<N0.A  rwo^orYTP  ' ^ 

|»A»G.TIX  ?.f4 
[ »ASG.T6mI  »,T,DAf3e 
I »ASG  TM  >A0C.T>»3W 
•RUM 


Figure  17.  Deck  structure  for  COPYTP. 


b.  Input  Data  Cards.  As  input,  COPYTP  requires  a data  deck  consisting  of  one  card 
for  each  segment  of  the  ephemeris  to  be  punched.  Each  data  card  has  the  format 
shown  in  table  14.  The  maximutn  number  of  data  cards  is  set  at  100.  The  data  cards 
can  specify  a mixture  of  satellites  and  ephemeral  times.  C(JPYTP  reads  the  first  record 
of  the  tape  and  determines  the  starting  epoch  and  satellite  on  the  tape.  The  epoch  and 
.satellite  of  each  data  card  is  compared  to  the  tape  parameters  and  the  cards  for 
satellites  and  times  not  on  the  tape  are  skipped.  The  remaining  cards  arc  sorted  by 
increasing  epoch.  The  requested  arcs  are  then  sequentially  extracted  from  the  tape. 


I 


CARD  NAME 
DATA  CARD 


Tabic  14.  COPYTP  Data  Card  Format 


PARAMETI-K 

CARD 

COLUMN 

FORMAT 

Termination  Symbol 

1 

11 

= 0 normally 

= 9 last  card  of  data  deck 

(Jrbit  Number 

2-7 

A6 

Satellite  Identification 

8 

A1 

Ephemeris  Segment  Epoch  (GMT) 

Month 

12-13 

12 

Day 

14-15 

12 

Year 

10-17 

12 

Hour 

18-19 

12 

Min  utc 

20-21 

12 
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20.  OUTPUT  ERROR  MESSAGES.  <i.  SATELLITE  XXX  NOT  FOUND  IN  TABLE 
(DATA  CARD)— The  satellite  code  in  the  pass  ID  is  in  error.  The  correct  satellite 
identification- satellite  number  correspondence  is  defined  as: 


Satflliio  Number 
Card  ('ode  on  T.rpe 

W 

X 2807 

V 2965 

b.  * ’ ' ERROR  ***  (DATA  CARD)— The  epoch  requested  on  the  data  card  is 
earlier  than  tlie  starting  epoch  on  the  tape. 

I NO  DATA-  The  data  has  been  sorted  and  none  of  the  requested  epheineris 
segments  are  on  the  tape  that  has  been  mounted. 

</  **•  PASS  XXXXXX  COULD  NOT  BE  LOCATED-Thc  requested  segment 
should  be  on  the  tape:  however,  it  could  not  be  found. 

e.  EOF  ON  TAPE-End  ot  Hie  encountered  on  tape— last  record  on  tape  is  printed 

i out. 

f.  INPUT  ERROR  XX  STATEMENT  Y— Type  XX  input  error  occurred  in  a record 
to  be  skipped  at  statement  Y in  the  program.  Recoverable. 

l>.  INPUT  CARD  ERROR-A  LAPUN  data  card  is  in  error— card  is  skipped. 

h.  CARD  EOR— The  termination  (end  ot  record)  symbol  is  missing  in  the  LAPUN 
data  deck. 


SECTION  V.  SUMMARY 

The  rXIPPLR  Point  Positioning  Program  h.is  been  in  production  use  at  the  DMA 
Topographic  Center  for  tlie  past  5 years.  It  has  proven  to  be  a reliable  software 
program  and  takes  advantage  of  all  the  inherent  accuracy  present  in  the  Doppler 
observations  from  the  portable  Doppler  receivers.  It  has  also  been  used  as  a simulation 
tool  to  determine  expected  geodetic  accuracies  under  various  operational 
configurations  and  constraints.  This  capability  has  proven  to  be  useful  in  planning 
Doppler  operations  to  support  both  mapping  and  precise  geodetic  applications. 
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